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Abstract. We discuss a microscopic scheme to compute the rigidity of glasses or the plateau modulus of supercooled liquids 
by twisting replicated liquids. We first summarize the method in the case of harmonic glasses with analytic potentials. Then 
we discuss how it can be extended to the case of repulsive contact systems : the hard sphere glass and related systems with 
repulsive contact potentials which enable the jamming transition at zero temperature. For the repulsive contact systems we 
find entropic rigidity which behaves similarly as the pressure in the low temperature limit: it is proportional to the temperature 
T and tends to diverge approaching the jamming density with increasing volume fraction as limj^o fl/T oc l/(0j — <p), 
which may account for experimental observations of rigidities of repulsive colloids and emulsions. 
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1. INTRODUCTION 

Supercooled liquids, glasses and jamming systems ex- 
hibit rich visco-elasticity [1]: at shorter time scales called 
as the j3 -regime such a system behaves as a solid with 
finite rigidity while it behaves as a liquid with high vis- 
cosity at longer time scales called as the a -regime. These 
features appear clearly in the relaxation of shear-stress 
(see Fig. 1) which follows after switching on a small 
shear-strain y (See Fig. 2). Approaching the glass tran- 
sition point the separation of the time scales between 
the two regimes become enormous so that a supercooled 
liquid behaves essentially as a quasi-static solid for a 
very long time. The shear-modulus or the rigidity is 
the most basic quasi-static property which distinguishes 
solids from other states of matters. 

Among the various types of glasses a class of systems 
like densely packed repulsive colloids, emulsions, foams 
and granular particles [2] exhibit an interesting common 
feature called as the jamming transition: the characters 
of the amorphous solid state change around the so called 
jamming point at a certain volume fraction at zero 
temperature. This is manifested in various quasi-static as 
well as certain dynamic properties of such amorphous 
solids [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
18, 19, 20,21,22, 23]. 

The wide time separation between the two regimes 
suggest that the quasi-static responses in the j3 -regime 
may be analyzed by some statistical mechanical ap- 
proaches such as the cloned liquid approach which com- 
bines the replica method and liquid theory [24, 25, 13, 
26]. The latter is a first principle, microscopic approach 
within the framework of the so called random first order 
transition (RFOT) theory [27, 28]. Indeed we showed re- 
cently that the replica method provides a trick to disen- 
tangle the j3 and a-like responses and extract the quasi- 
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FIGURE 1. Schematic picture of stress relaxation in super- 
cooled liquids. A simplest protocol is considered: switch on 
the shear-strain of amplitude y (See Fig. 2) at time t = and 
measure the relaxation of shear-stress o"(f ) which follows. Here 
we define fl(t) = a(t) jy. The plateau value G p of fl(t) is called 
as the plateau modulus which represents the effective rigidity or 
shear-modulus of metastable glassy states. 



static part of the responses in supercooled liquids and 
structural glasses [29, 30]. 

The purpose of the the present paper is two fold: we 
first review the method [29, 30] developed for systems 
with analytic potentials such as the Lennard-Jones poten- 
tial which are suited for usual molecular glasses. These 
systems behave as harmonic solids at low temperatures, 
i. e. systems of random spring networks. Then we at- 
tempt to extend the method to account for systems which 
exhibit the jamming transition. The essential ingredient 
in such systems is the repulsive contact potentials such as 
hard-spheres and some soft repulsive contact potentials 
which are not analytic. Based on this method we analyze 
the behaviour of the rigidity of hard-spheres and soft re- 
pulsive contact systems in the low temperature limit ap- 
proaching the jamming density 0j from below. It appears 



that our result accounts for some experimental observa- 
tions of the rigidity of repulsive colloids and emulsions 
[4, 5, 6]. 

The organization of the paper is as follows. In the next 
section we introduce the two distinct classes of systems: 
systems with the analytic potentials and the repulsive 
contact potentials. In the subsequent sections we review 
our strategy [29, 30] to extract the quasi-static response 
functions of supercooled liquids and glasses based on the 
cloned liquid approach. Then we review the basic fluctu- 
ation formulae of the rigidity and our previous scheme 
to compute the rigidity of harmonic glassy systems re- 
ported in [29, 30]. Finally we discuss the extension of 
the method to the cases of repulsive contact potentials 
and analyze how the jamming transition is reflected on 
the entropic rigidity approaching <j)j . 



2. MODELS 

We consider a generic system of N particles (z = 
1,2,..., AO in the 3 -dimensional space with volume V 
interacting with each other through a two-body poten- 
tial v(r) which only depends on the relative distance r 
between particles. The potential part of the Hamiltonian 



can be written as, 



U : 



Iv(ry) 

(H) 
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where (ij) stands for summation over the N(N - l)/2 
pairs of the particles and r, (i= 1,2,..., N) represents the 
position of the particles. We suppose that the temperature 
T is low enough and the number density p = N/V is high 
enough such that the system is in a supercooled liquid or 
a glassy metastable state. 

We consider two distinct classes of systems: 

• Harmonic systems: the potential v(r) is an analytic 
function of r for r > like the Lennard-Jones po- 
tential. Presumably this class of systems is relevant 
for molecular glasses. 

For an explicit model computation we consider the 
soft-core potential v(r) = e(r/a)~ 12 where e and a 
are the unit energy and length respectively. 

• Repulsive contact systems : the potential v(r) is re- 
pulsive and has a definite cut-off at the scale the par- 
ticle size a like the hard-spheres. Presumably this 
class of systems is relevant for emulsions and repul- 
sive colloids. The density is a crucial parameter in 
these systems and it is convenient to represent it via 
volume fraction <p which is related to the number 
density p as = (n/6)a 3 p. An important feature 
of this class of systems is that they exhibit the jam- 
ming transition at T = by increasing the volume 
fraction <j) up to some jamming density 0j. 



For an explicit model computation we consider the 
soft-particle potential v{r) = e(l — r/a) 2 9(l —r/a) 
where 9 (r) is the step function. In the present paper 
we limit ourselves to the range of volume fractions 
(j) < where the system behaves as hard-spheres in 
the zero temperature limit T — > 0. 



3. DISENTANGLEMENT OF THE 
INTRA STATE AND INTER-STATE 
RESPONSES 

3.1. A mean-field picture: ensemble of 
metastable states 

Let us take the basic energy landscape picture [31, 32] 
: we consider that the equilibrium state of a supercooled 
liquid or a glass can be described in terms of a statis- 
tical ensemble of metastable states, which might be in- 
terpreted as metabasins [32] each of which is a union of 
inherent structures [31]. Let us label the metastable states 
as a = 1 , 2, 3, . . . and denote the free-energy per particle 
of the a-th state as f a . Then the equilibrium free-energy 
of the system may be expressed formally as, 



F(h) = -k B TlogZ(h) 
Z(h)^e 



-Wa(h) 



(2) 



where k B is the Boltzmann's constant and j3 = l/k B T. 

We have introduced a parameter h which represents 
a generic infinitesimal probing field, such as the shear 
which we will focus on in the present paper. The linear 
susceptibility to the external field h can be seen to take 
the following generic form, 



1 d 2 F(h) 



N dh 2 
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with 



X=llXa}} X = mi[o 2 J-[[o a f) (4) 



Here o a = df a (h)/dh\ h=0 is the equilibrium value of an 
observable o which is conjugated to the external field h. 
Similarly % a = do a /dh = d 2 f a (h)/dh 2 \ h=0 is the as- 
sociated linear susceptibility within a given metastable 
state a. In (4) [[. . .]] stands for averaging over the ensem- 
ble of the metastable states defined as, 



£ gg -^/°(°)... 
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The important feature evident in (3) is that the total lin- 
ear susceptibility % is the sum of two distinct parts: % and 



X associated with the response within metastable states 
and response due to jumps between different metastable 
states. Physically % can be regarded as the quasi-static 
response within the j3 -regime and % can be related to the 
response in the a-regime. Although the two parts have 
very different characters, (3) implies they are mixed up 
in the total response. We wish to disentangle the two. Let 
us discuss below how the replica trick works for this pur- 
pose. 



3.2. Response of a cloned system 

Let us consider a cloned system [33, 34, 24] which 
consists of m replicas of the same system labeled as 
a = 1,2, ... ,m. The free-energy of the cloned system is 
defined as, 



F m ({h a }) = - — \°gZ m ({ha}) 

m 

Z m ({h a })^e- N ^M" a ) 



(6) 



Note that there is only one summation over the 
metastable states instead of m summations. This means 
that we are assuming that m replicas are not allowed to 
fluctuate independently from each other but forced to 
fluctuate together over different metastable states. Yet 
the replicas are allowed to fluctuate differently from 
each other within the metastable states. How to realize 
such a situation in practice is a non-trivial task by itself 
[33, 34, 24, 13] as we discuss shortly later. 

The key point is that we have put different probing 
fields h a (a = 1 , 2, . . . , m) on different replicas in (6) [29, 
30]. It naturally lead us to define a sort of generalized 
linear-susceptibility of a matrix form, 
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_ 1 d 2 F m ({h}) 



N dh a dh h 
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= Xm8ab + Xm- (7) 



where % m and % m are almost the same as % and % defined 
in (4) but evaluated by replacing [[. . .]] by [[. . .]],„ defined 
as, 



Z,„(0) 



(8) 



Quite remarkably the 2nd equation of (7), which can be 
easily verified, implies that the j3 and a-like responses 
can be distinguished from each other: 



% = lim Xm 



lim X 



Xm ~ Xaa Xa^b 
•n Xm — Xa^b 



(9) 
(10) 



3.3. Cloned liquid 

Here let us briefly sketch how to implement a cloned 
system [24, 13]. The basic idea is to introduce a system 
of an artificial molecular liquid in which each 'molecule' 
i = 1,2,...,N consists of m particles belonging to differ- 
ent replicas a = 1,2, ... ,m. The particles are allowed to 
fluctuate only within the molecule of size A, which is in- 
terpreted physically as the cage size. The cage size A is 
determined by a variational principle (see below). Exis- 
tence of a solution with a finite cage size A < °° implies 
existence of metastable states [34, 24]. 

The coordinates of the particles rf (i = 1 ,2, . . . ,N) can 
be decomposed formally as, 



rf=R,+< R, 



1 



(11) 
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where uf stands for fluctuation of the particle belonging 
to the fl-th replica with respect to the center of mass R, of 
the molecule. The fluctuations within the molecules are 
assumed to obey the Gaussian statistics with the mean 
and variance given by [24], 

(«fW = ((utr(u b j r) cage =A(8 ab -l i )8^8 i j 

(12) 

Here fi (and v) represents a component of 3-dimensional 
vectors jl = x,y,z. The factor 8 a b — ^ reflects the con- 
straint £™ =1 (u a ) = 0Z. 

The free-energy of the molecular or cloned liquid 
G m {A) of a given number of replicas m and the cage size 
A can be obtained as follows [24]. First one integrate out 
the fluctuations within the molecules which amounts to 
replace the original interaction potential v(r) by a remor- 
nalized one v e s(r,A). [24, 13] In the case of analytic po- 
tentials it reads as [24], 



Veff(r) = v(r) - (1 - m)-jV 2 v(r) + ■■■■ 



m- 



(13) 



Then one is left to integrate out the CM positions of the 
molecules interacting with each other via v e ff(r,A) and 
subjected to a heat-bath at an effective temperature T/m. 
Eventually we have to take the m — > 1 limit (see (10)). 

The strategy is to start from sufficiently small m(<IC 
1) so that the cloned system remain in the liquid state 
because the effective temperature T/m becomes suffi- 
ciently high even if the actual temperature T itself is very 
low. Then standard density functional methods of the liq- 
uid theory [35] allows one to compute the free-energy 
G m (A). The value of A is determined by minimizing the 
variational free-energy G m (A) with respect to A yielding 
F m = mi\\AG m {A). Let us denote the value of the cage 
size A at the minimum as A* (m, T,p). 

The last step is to take the limit m -> 1~. It turns out 
that at temperatures T below the ideal glass transition 



temperature, i. e. the Kauzmann temperature 7k (p), one 
finds a characteristic value m* — m*(T,p) in the range 
0<m*(r,p) < 1 such that 



lim F m (T,p)=F m *(T,p). 

m-5- 1 - 



(14) 



holds. The Kauzmann transition temperature 7k (p) can 
be obtained by solving m* (7k (p),p) = 1- The reason 
behind (14) is actually the entropy crisis mechanism, i. e. 
the ideal glass transition, taking place along the m-axis. 
We refer the readers to Refs [24] for the details. The 
above observation implies that j3 and a-like responses 
(see (10)) can be obtained as, 



X — X" 



X — Xii 



(15) 



4. MICROSCOPIC COMPUTATION OF 
THE RIGIDITY OF GLASSES 

4.1. Static linear response to shear 

Now let us study static linear response to shear. For 
clarity we recall and compare the well known fluctuation 
formulae for the shear-modulus for analytic potentials 
[36] and hard-spheres [37] in the case of simple shear. To 
this end we consider a system of particles i = l,2,...,N 
put in a rectangular container of volume V and perturb 
the system by a shear-strain of infinitesimal amplitude y 
on the container (See Fig. 2). Most important feature of 
the shear is that it just changes the shape of the container 
but not the volume V (and thus the density p). 

The free-energy F(y) of the system may be formally 
expanded in power series of y as, 



F(y)/V = F(0)/V+yo + ^h- 



(16) 



Here the coefficients of the 1st and 2nd order terms in 
the expansion defines the shear-stress a and the rigidity 
or the shear-modulus ji. Microscopic expressions of a 
and ji can be obtained as follows. 

The free-energy F(y) can be expressed formally as, 



F(y) = -k B ThgZ(y) 
Jr( 7 ) i=l A th 

N j3j 

ff±Jl e -PL { ij)v{rij) 



the 2nd equation we changed the integration variables 
from r to r' (See Fig. 2) which allows us to change 
the integration region y(y) back to the unperturbed one 
^(O). The 2nd equation allows us to easily obtain the 
expansion of the free-energy F(y) in power series of y. 

Z 



Here A t h is the thermal de Brogile length. The subscripts 
y{y) represent the range of integrations, including not 
only the volume V (which is invariant under shear) but 
also its shape, parametrized by the shear-strain y. In 
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FIGURE 2. The geometry of a sheared container. 

Shear-stress. Microscopic expression of the shear- 
stress is obtained as, 



1 dF{y) 



V dy 



where 



7=0 ' (ij) 
is the thermal average, 
1 



I>(r y )) (18) 



1 1 1 3 



= 1 \ 



y "' z(o)J m i 

Here we introduced the local shear-stress. 

c(r) ee zxrv'(r) 



(19) 



(20) 



with short-hand notations x = x/r, y = y/r , z = zjr 
where r = y^x 2 +y 2 + z 2 - 

Remark: note that the expression (20) is ill-defined 
for hard-spheres. However the average shear-stress 
Eq. (18) can still be evaluated safely by noting that 
-pV '(r)e~0 v M = (e~0 v M)' becomes a delta function 
8(r — a) for hard-spheres and soft repulsive contact 
systems in the T —> limit. 



Shear-modulus. Similarly the microscopic expres- 
sion or the static fluctuation formula of the shear- 
modulus is obtained as, 



(17) 



1 d 2 F(y 



V dy 2 



(21) 



y=0 



-P>E E «cr(r w )(T(r m „)) - (a(r kl )){a(r mn ))) 

(kl) (mn) 



where ji\, om is the so called Born term defined as, 
1 



V 



£f 2 [r 2 v"(r)i 2 + /-v'(/-)(l-i 2 



(22) 



The Born term jllborn represents the instantaneous, affine 
response to shear while the 2nd term in the r.h.s of the 
2nd equation (21) represents the so called non-affine 
correction due to stress relaxation [36, 30]. 

The above expressions (21) (22) are problematic for 
the repulsive contact systems. Especially the Born term 
is formally infinite for hard-spheres. The Born term 
Eq. (22) stems from direct spatial derivatives of the lo- 
cal stress (20) which does not exist in this class of sys- 
tems in sharp contrast to the harmonic systems (See the 
remark below (20)). Then for this class of systems it is 
more convenient to use an alternative but equivalent ex- 
pression [37], 



as well as the interaction potential as, 

v(if y )=v(J?y)+Vv(r)U^.(u? 
d 2 v(r) 



+ ^E 



dr^dr 



Km.... 



r=R h 



As the result of the cloning we obtain the rigidity 
matrix in the form of the anticipated matrix structure (7), 



Hab = A<U + £ 



(24) 



with the rigidity of metastable states (plateau modulus 
G p ) obtained up to 0(A) as, 



-(H) 

-£ £ «j3a(r H )j3a(r B 

{kl) (mn)^(kl) 



(23) 



eft 
reborn 



(25) 



t ))-(Po(v kl ))(l3o(r mn )))} 



which can also be obtained from (17). In the derivation 
one has to perform some integrations by parts in order 
to get rid of the Born term. In (23) we dropped off 
some terms which cancel out with each other exactly 
in isotropic systems. This assumption is valid for the 
systems we consider below. 

In liquids the terms on the r.h.s of (21)(23) cancel with 
each other to realize jU = which reflects the transla- 
tional invariance of liquids. 



5. HARMONIC SYSTEMS 

Let us now discuss the rigidity of the glassy states of the 
harmonic systems whose potentials v(r) are analytic for 
r > 0. Here we present a summary of the results reported 
[29, 30] which will be compared with the case of the 
repulsive contact systems analyzed in the next section. 
First we consider a 'free' m-replica system without any 
'cloning': the m replicas are totally independent from 
each other. For such a system we can naturally define 
a rigidity matrix p, a i, and find its microscopic expression 
similarly as to the one for the single system (21), 

d 2 F m 



P-ab 



I 



V dy a y b 



Pbom&ab 



-0VE E ((a(r a kl )a(r b mn )) (a(v a kl )){a(v b mn ))), 

(kl) (mn) 

where jUbom and the local shear-stress <j(r) are the same 
as those defined in Eq. (22) and Eq. (20). 

Now we switch on the cloning following the prescrip- 
tion discussed previously using the decomposition of 
the coordinates (11). This allows us to expand the local 
shear-stress o(rfj) as, 

a(rf j ) = a(R ij ) + Va(r)\ r=R .(u?-„») + ... 



m* 



P Jd 3 r\(Pm*Va(r))\ 2 g T/m , !p (r) 



3 3 

a r\d r2 



(/3m*Va(n)) • (j3m*Va(r 2 )) (g 3 )r/™.(ri,r 2 ) 
while the a-like part of the response is given by 



M = --r 



m* 



(26) 



because of the sum rule Y!b=i V-ab — reflecting the plain 
fact the cloned liquid as a whole is just a liquid. Phys- 
ically the latter suggests static analogue of the yield- 
ing processes (See [30] for discussions). In the above 
equations m* = m*(T,p) and A* = A*(m*(T,p),T,p) 
are the values determined in the course of the evaluation 
of the free-energy of the cloned liquid discussed before. 
Here gr,p(r) and (g3)r,p(r 1,1*2) are the radial distribu- 
tion function and the three-point correlation function of 
the liquid at temperature T and density p respectively. 

The term jJ.^ m represents the Born term (22) associ- 
ated with the renormalized potential (13), which itself 
consists of the original born term (22) at O(A ) and cor- 
rections at 0(A). The last two terms on the r.h.s of (25) 
represents the effect of the stress relaxation due to the 
fluctuation inside the cages, i.e. j3 -relaxation (See Fig. 3). 

We show the result of a model computation on the 
binary soft-core system v(r) <~ (r/a)~ 12 in Fig. 4. The 
result at finite temperature T > was already reported 
in [29, 30]. Here we added the result of a computation 
performed directly in the zero temperature limit T — >• 0. 

The cloned liquid computation can be performed in 
the zero temperature limit T — >• by introducing scaled 
variables a ~A/T and T — T /m [14]. The expression of 
the rigidity jl given Eq. (25) also suggest it has a limiting 
value in the T — > limit. 

For the binary soft-core system we consider T -> 
limit of the the theory formulated in Ref [25] which 




FIGURE 3. Schematic picture of stress relaxation due to 
motions of particles inside cages. The local motion of the 
particle i can relax local stress between the particle i and the 
surrounding particles j\,j2, 



employs the binary HNC approximation for the com- 
putation of the liquid free energy and the radial dis- 
tribution function gr(r). As the result we find the ef- 
fective temperature of the cloned liquid converges to 
X = limj^o T/m*(T) ~ 0.115 within the 1st order cage 
expansion. This value is close to the Kauzmann transi- 
tion temperature 7k /£ ~ 0.12 where m*(7k) = 1 [25]. 

The result shown in Fig. 4 suggest the rigidity at low 
temperatures T < 7k behaves essentially as if at T = 0. 
Moreover the numerical values of the plateau modulus 
G p observed by MD simulations [38, 39], and the value 
of the rigidity of the inherent structures [40] compare 
well with the theoretical values shown in Fig. 4. These 
observations strongly support the usual view that the 
metastable glassy states at low enough temperatures are 
harmonic solids which can be efficiently described as 
systems of random spring network (see also [41] for a 
related work). 

Note also that the rigidity becomes significantly 
smaller with increasing temperature above 7k. We in- 
terpret this as reduction of the rigidity [30] due to ther- 
mally activated plastic events among a union of inher- 
ent structures [31] belongin to a common metabasin [32] 
or a metastable state. The rigidity apparently vanishes 
(crosses 0) around T/e ~ 0.22, suggesting melting of 
metastable states, which happens to be rather close to the 
so called MCT critical temperature T c /e ~ 0.19 - 0.22 
[42] [26] . However the first order cage expansion which 
we have employed does not allow us to locate the MCT 
transition temperature at which the glassy solution with 
finite cage size A < °° disappear presumably by a spin- 
odal like mechanism [34]. As we noted in [30] we rather 
consider at the moment that, at least in the mean-field 
sense, the rigidity should exhibit a discontinuous be- 
haviour. This is because the expression (25) actually 
implies that the rigidity is a function of the the cage 
size A which is predicted to behave as A(T) — A(T C ) 



\/T c — T approaching the dynamical temperature T c from 
below by the mode coupling theory (MCT) [43]. Indeed 
such a discontinuous behaviour of the rigidity has been 
suggested in an alternative formulation of the replica 
approach [44]. 
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FIGURE 4. The rigidity of a binary soft-core system. The 
born term which represents the instantaneous, affine response 
and the rigidity jl which includes non-affine corrections by 
stress relaxation due to fluctuations inside cages are shown. 
The results are T > are reported in [29, 30]. The Kauzmann 
temperature is Tk/e ~ 0.12. 



6. REPULSIVE CONTACT SYSTEMS 

Finally we are in the position to analyze the rigidity of 
the glassy states of repulsive contact systems: the hard- 
spheres and generic soft repulsive contact systems in the 
low temperature limit T — >• 0. 

Again we first consider the rigidity matrix of the 'free' 
m-replica system, which is obtained as, 



£<j8a(4))<j8a 6 (4)> 

l{kl) 



£ £ ((Mr^/MOM/MrDXMd,))) 

(Id) {mn)^(kl) 

As expected it is similar to the expression Eq. (23) for 
the single system. 

Then by switching on the cloning we obtain the rigid- 
ity of metastable states (plateau modulus G p ) up to 0(A) 
as, 



j3£ 



1 /A* 

m* \m* J n a 3 
(Vj8m*a(r 1 ))-(Vj8m*a(r 2 ))(g 3 )r/m*,0(i-i ! i-2) (27) 



- / d 3 nd 3 r 2 



Most important difference from the case of the harmonic 
systems (25) is that the Born term is apparently absent 
here. Consequently the rigidity is at order 0(A°) and 
starts only at 0(A). The relevant term at 0(A) is again 
the one related to the thermal fluctuation of the shear- 
stress due the the fluctuations inside cages (See Fig. 3). 

In order to make further progresses, we approximate 
the three-point correlation function g3(ri,r 2 ) by the 
Kirkwood approximation, 

(g3)r,4>(ri,r 2 ) ^ gT^(n)gT^(r 2 )g T ^(\Y]_ -r 2 |). (28) 

Then using the cavity function yr,^(r) defined as 
grj (r) = y T 4 (r)e - P v M we find, 



nf \m I % O 



(29) 



with 



„HS 



c~-j d\,d\ 2 yT/ m *A\^-^\)e-^ v ^- r ^ 

( de-P m * v ^)\ 

/ de- pm * v ^ \ 

• v 2 \ yTim*$(n)z2hn — — — I 

-(yf(a)) 2 J 'ddidCh jdnd^rjri 

3^(|ri-ri|)0(|n-ri|-fl) 
Vi (zixin8(n - a)) ■ V 2 (z 2 x 2 r 2 8(r 2 - a)) (30) 

In the last equation we took T — > limit which greatly 
simplifies the calculation. The evaluations of the inte- 
grals over the polar coordinates are tedious but straight 
forward. First the integrations along the radial coordi- 
nates n and r 2 can be done exactly via integrations by 
parts. Subsequently the integrations over the solid angles 
Cli and £l 2 can be simplified, with the help of spheri- 
cal harmonics, to simple one dimensional integrals over 
x = cos(0i 2 ) where 9\ 2 is the relative angle between the 
two solid angles. If we make a further approximation 
3^ S (l r i — r 2|) — y$ S ( a ) the last integrations can also be 
done exactly and we finally obtain, 



(31) 



The above results suggest finite rigidity jl of the repul- 
sive contact systems which is proportional to the temper- 
ature T meaning that it is of entropic origin in sharp con- 
trast to the harmonic systems discussed previously whose 
rigidity is essentially mechanical. 

In order to study how the rigidity depends on the vol- 
ume fraction 0, we need to know the values of m*(T = 



0,(j)),A*(T = 0,0) andy" s (a). Fortunately they are pro- 
vided by recent studies on the hard-sphere glass [13] 
and a soft repulsive contact potential system [14] close 
to the jamming density, more precisely the so called 
glass close packing density 0gcp- According to the latter 
works m*(r = 0,0) ~ci(0 GC p-0) andA*(r =0,0) ~ 
c 2 (0gcp — 0) approaching 0gcp from below with c\ and 
c 2 being some positive constants. The value of y^ cp ( a ) 
is also positive. Using these information in our result we 
obtain, 

lim pp. = "7~~~~~r (32) 
r-><r^ 0gcp-0 

with the numerical pre-factor given by, 



, 0GCP c 2 113 HS , ..3 
c = 6n— — T — (y » (a)) . 



(7 3 c\ 120 v/fccpV 



(33) 



Using the numerical values of constants reported in [14] 
we find c ~ 0.7. A remarkable feature is that the scaling 
of the rigidity jl found above Eq. (32) is exactly the same 
as that of the pressure p found in [13] and [14]. We note 
that our result is different from that of Ref [10] which 
predicts somewhat stronger rigidity jl « p 3 / 2 based on an 
effective (yet microscopic) harmonic description devel- 
oped for the inherent structures of the hardsphere glass. 

The scaling (32) agrees with the plateau modulus G p 
observed by MD simulations of the stress relaxation on 
the same system [39, 45]. However the evaluation of the 
numerical factor c (33) may be improved in several re- 



spects: 1) the approximation 



„HS/ 



ri 



„HS 



(a) can 



be avoided by doing numerical integrations 2) higher or- 
der corrections terms of the cage expansion due to renor- 
malization of the potential [13, 14] can be considered 
3) better evaluation of the 3-point correlation function 
<?3 ( r 1 j r 2 ) than the Kirkwood approximation may be con- 
sidered. We have checked that the item 2) amounts to 
reduction of the value of c by an amount of about 20% 
which will be reported elsewhere. The item 3) would be 
challenging but worthwhile. 

It is interesting to compare the above result with 
some experimental observations of the rigidity of densely 
packed repulsive colloids [4] and emulsions [5, 6]. The 
experiments were performed at the room temperature 
which is actually a very low temperature for these sys- 
tems. For example the reduced temperature can be esti- 
mated as k B T/e <~ 10~ 5 for the emulsions system [5, 6]. 
However the experimental data reveal presence of finite 
entropic rigidity which rapidly increase approaching the 
jamming density from below. A striking feature found by 
the experiment on the emulsion system is that simultane- 
ous measurement of the pressure reveals that the pressure 
and the shear-modulus behave very similarly (See Fig. 3 
of [5]). Thus it appears that the our theoretical result is 
consistent with the experiment. 

In the present paper we have limited ourselves to the 
volume fractions < 0gcp> but it is straightforward to 



extended the present approach to the jammed region 
> 0gcp concerning the systems of the soft repulsive 
contact potentials. It is important to note that for this 
class of systems one cannot rely on the usual picture 
of harmonic solids naively [15, 16, 17]. We will report 
the results elsewhere with detailed comparisons with the 
known results [3, 5, 6, 7, 8, 10, 12]. 

7. CONCLUSIONS 

In the present paper we first reviewed a microscopic ap- 
proach to study the rigidity of structural glasses or the 
plateau modulus G p of supercooled liquids based on the 
cloned liquid approach. Then we discussed how to ex- 
tend the method, which has been limited to the cases of 
harmonic glasses, i. e. systems with analytic potentials, 
to the cases of the glassy repulsive contact systems like 
the hardsphere glasses or soft repulsive contact systems 
in the low temperature limit. We found the entropic rigid- 
ity of this class of systems exhibit divergent behaviour 
much as the pressure approaching the jamming density 
from below, which appear to be consistent with experi- 
mental observations on repulsive colloids and emulsions. 
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